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Abstract 

For llnew wave prtp^gatjon io anisotropic media, the pfindpk of 
superposition is stilt hoid. Ttit deoonposiUoo of ^mbtV dipole 
log is based on tbifi priocipleL In (he forward decompttition in- 
lioe and erossiine acoustic dipole ioga u aoy ayimutbal aogte b th« 
proJectioQ of mcasuremefttfi along the priociple direction Of tbe focw 
matlot). Id tho invee&e dceomposiUan the measureiBeqtt along the 
principle direction can be cOAlLructed from the Ofthogonai pair of 
inline aad croaAline acou$tjc dipole log. The analytic formulas for 
IxHh forward and iovecM decojn position of dipole log arc derived 
in this paper. The inveiae decompodtion formula Is the «Otutioa 
in the leaatro^uare seiwe. Numeficai examples are demodstraied 
for acoustic dipole log decompositioa io isotn>pie and vii«otrOpM: 
formitions. The synthetic dipoie log is calculated by 3D finite 
difference method. The namericaJ examples (ir« abo aJsowo that 
the inveise decomposition formula nfOrks very well with n^y data. 
This inverte decomposition formula will very uaerul to procesS the 
field acoatilic toeing daU in anisotropic formation. It caa provide 
the direction of the formation anisotropy well aa Che degr^ of 
anisoiFopy. Becaiae acoustic dipole logging is in ths near iield dis- 
tance, the partWc motion id very complicated. Only In tbe principle 
direction the partjde motion b linear polarized. The Initial parti- 
cle motion with dipole dOUrce at arbitrary azimutha] aAgle tends 
to point In the fast shear wave directioo. But it wiD be difficult 
to use this Information to ^d »tabiB esUnation of fast shear wave 
direction. 

Introduction 

The crust of tho Earth ta slightly anlsotropk, wiikh la rdated (o 
geological processed. For example, antaotropy can be caused by 
aCgned fraauros in the rock. Pine layered sedimentary cocks also 
poeBBss tranaverse isotropy. AoDustic logging provides a tcchniqae 
to meaiara ^mation anbotropy froa drilled weU.' 

To gnderstand the effects of aabotropy on the fluid-Qted bore- 
hole wave propagation ace critical to proccM acoustic logj^ng data 
in anisotropic IbrmattOA. Several numericaJ techniques have been 
developed to simulate ■couatie logging, $«ch as discrete wavennm- 
ber method (While aod Tbagtaow, Idgl) Chan and Tsaagt 1983; 
Schmitt« 1989), pertarbatloD method (Ellefsen, 1990; Sinha et al., 
19d4) and anit« difference method (Leslia and AandaU. 199% Chen 
et al., 1995). Fbr iranavei^ isotropic formaiioa with tbe symme* 
try axb aligned vuHh the borehole axis, the solttUon b wdl known 
and the different type of waves are well understood. But when the 
aymmetcy ans of the TI fbrmation makot aq vbitrary angle with- 
the borehole akb, tbe borehole wave propagation proUem becomes 
much more compUcaied. Most of our knowledge are come from the 
perturbation solution. On the other hand field data obsenAtions 
in the anisotropic formation are available, such as borehote waves 
generated from vertical seismic profile (Barton aod Zoback« 1988; 
LeveQle and Seriff, 1989) and borehole dipole logging (Eameraoy, et 
al., 1994). 

Th« purpoca «i thU r»p«r u t#y ti> »n«w«r nk* Ia1L(»wIa|^ 
quesUona. How to construct the inline and cros&llne dipoie log sU 
any aatimuthai angle from the measurement abng the principle axb? 



What the particle motion of 4ipole log k)ok& like in an anisotropic 
(brmadon? Kow CO raeover the measurements along the principle 
axb from the inline and croBsline logs at arbitrary actmuthal an^. 
The first and the third question are the forward aod tbe inverae 
problem pair. In order to do numericja experimenta the 3-D time 
domain finite difference method b used to generate synthetic dipole 
bg. The formation b assumed tranaverwty isotropic «ath ita aym- 
metry axis perpendicular to the borehole axb. 

Ebrwvrd Deeompositioh of Dipol« Log 

Although we consider wave propagation in wiuotropk media, the 
wave motion b vtiU governed by linear differential equations. Tbe 
principle of 6up«rp««tlon b valid. We assume that the X and the 

Y axb 4re aliened vrath two prindple axes of anisotropy. The dipole 
Gource b In the X-Y plane and Gij b the Oreeft*a function (index 
I and i reprettAt x,y). The general form of dipole meaauiement is 
^^Gii* where R represent receiver and % repreacnts source. 
The symmetry of the model reqnires a G^ ^ i^. We de- 
fine that s Gtt • S b the inline dipole measareroent In the X 
direction and Oy = * inline dipole measurement in the 

Y direction. These measurtmente can bu computed by Rnite dlf' 
(erenee method. The same dipole source \;^th azimuthal angle 9 
(counterclockwise from the X axb)» can be decompOted into the X 
and the Y directions. Then the mcasucement along the X axb will 
be: 



and the measuremeni along the V axb wiB be: 



(2) 

' fa) order to demonstrate we use 3D finite difference method to 
compute synthetic dipole log in isotropk: u wdl as anisotropic 
formations. It is a fourth-order scheme and implemented on a 
massively parallel computer. A more detailed descripOons of the 
method can be found in the paper by Cheng et al., (1995). We 
'fim consider botropie formation, it b a slow formation with ^ 
wave velodty 9770 m/s. 5 wave velodty 1229 m/s. The water Btled 
borehole is 10 cm in radiss. The source b Kelly wavdet at cen- 
ter, frequency 2 kHs. The source and the receiver are located at 
tbe center of the borehote. The source and receiver separation is 
2 meters. These par&meteru are kept the same in the paper. The 
comparison of the direct finite difference computation with tbe pro* 
jection of Equation (1) b plotted in Figure II The dipole source b 
at a3Jmetbai angle o&45 degree.- The projected and the computed 
waveforms are in excellent agreement. 

• Two types of the T! formation are considered: (I) Bakfcen sh*4e. 
a fast formation (2) Austin chalky a slow formation. Tbe symmetry 
axb of the T] formation b perpendicular to the borehole axb. The 
elastic constants of these formations b give by Sinha et al. (1094). 
For the dipole source at aximuthal angel of 45 degree, the com- 
parisoos of direct computed waveforms with the projected one aro 

formation used is Bakken shale. Once again they are in excellent 
agreement. Equation (I) and (2) allow you decompose the dipole 
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Decomposition of dipoto log 



source into tho principle direct jofui kt\4 only compute the rappi»e 
at the principle diroctionA. U c&n save compuliog time for icous- 
tic dipole log sirnolatiort in an anisotropic formation with diffmot 
a^imuthd dlrectiDns. 

From two inline dipole raeasufemenla aJoag the X and the Y 
axis, the inline and £f03Sline dipole measurements ai aslmuthal 
angle i can bv eafiy obtained tfom the projection. We have: 



(3) 



vtrM4iint - -u*sin ^ + wjcea^ (4) 
fittb^tUute equation (l) and (2) uitoihe above eqoaliona waobtaSn 

*^ni<u4 = Cos' ^ + sin^ B 



and 



(5) 



= + Vf)slnB<u»$ (6) 
the S4me e<iu4tioDs a.re also given by Siph^ el al. ( 1904) and Esmec 
aoy 01 a|. (Id04), These utt our forward 4eeompo«tion formulaa 
for aquatic dipole log. To ehow some examples, projected the 
inUna and tho cKaanne dipok log in Austin chftJk formaUon ai dif- 
ferajit asimutbal angles. The waveforms ue plotted in Figure 2 In 
isotropic roripftUon the crosftUne dipole log is zero. In anisotropic 
formation th« crMollne dipole va<i(;&hea only when the dipole source 
IS aligned with tho symmetry axis. The croesline dipole log gener- 
ally is not zero, la logging data processing the crosHne dipole log 
CM be used todetermioe whether the formation shear wave velodty 
» anisotropic and the direction of prUcipIo axis of the ajiiaotr^py 
IR tlUIne dipole log (Figure 2) It Is dear from the waveforms that 
at SO degT«e asirouthai angle the arrival times are earlier than at 0 
degree. This \n the reflection of the shear wave velocity wimuthai 
dependent. The iniine dipole log can provide the measurements 
of anisotropy froqi the arrival time differences l>etweeQ dieenmt 
azimuthal — 



ParticU Motion &nd Polarl^tton Analysis 

The she^ wave splitting is often obiarved from the tehaeismic 
events. The parMde motion of shear wave is polarized along the 
principle direction. If juuHisiic dtpol« log \n anisotropic formation 
h also polarised along the principle directions, then we can find 
the principle direction from polariaaiion analysis. This will make 
it very eaay toaolve the inverse decomposition quotioQ. 

We define the mean value of N ahservationa of vanablo zii) as 
(Kanaaewich, IMl) * ' 

ID I 

and cawiaAM bctwnn two wisbUsT, »od x, » 



<?««(»„ „| = 1 , fi(„))((„((, _ E(,t)) (8) 

For the twocomponenu of dipofc log „,). ihecovariance matrix 
C can bo formed aa 

^' ^ cigenvalufls of the covariance matrix C and 



At > A), then the quantity 



P=l-:i2 



(10) 



estimate the rectillnearliy of the partide motion. When the rec- 
Ullnearity U high (A| > A3) the particic motion is close to linear 
polarisation. When the rectilineariiy ba near icro (A, :sf Aj) the par- 
ticle motion is near circular. The direction of polarizatiOQ may be 
meaeiired by considering the eigenvector associated with tiie largest 
eigenvalae. 

In Tl formations shear wave velocity is azimathaj dependent. 
The panicle moiions are plotted irt Figure 3 for Aoatin chalk. FVom 
the pardcle motion plots, the Uneir polarisalioos at a^imuthaj an- 
gles of 0 and 90 degree axe easy to understand. The jymmetry of 
the model requires that there is no motion perpendicular 10 the 
plane of the symmetry (G„ = Gy. = 0). The Hnear polarised par- 
tide motions at azimuchal angles <rf 0 and 90 degree also give some 
indirect aficutacy measure of finite difference method. When the 
dipole source is not aligned with the principle axis of the formation 
we would intuitively expect that it will polariacd along the two prin- 
dpie axes. But the plot of the particle motion did not show these 
polarisatioms dearly. The reason ts that oqr intuition are based on 
tho far-Qeld shear body waves in anisotropic media, Tlie dipole log 
Is using flexural wave which is a guided mode and the acoustic log- 
ging distance is belong to the near fieid. In general for one period 
of separation in time dom^n take? the distaoce about wavelength 
divided by the percent of anisotropy. For ezaaple a S% aoisotropy 
calces about 20 wavelengthes to gel one period separation id time. 

In order to verify the distaoce eiHsct we increase the source and 
receiver distance to 4.S meters In Austin chalk formation to do the 
finite dUttrence calcolaUons. Other parameters are not chasgod. 
The partide motion is plotted In Figure 4. It is very dear fcom 
the plot that the parUde niotkjD U started in the fast &hear wavt 
direction. Thb example supports that no clear polarized particle 
notion in dipole log b due to the nature of near-field loggiog. 

TherectjUnearity analysis ofthe d I pote particle motion in Bakken 
shale and Austin chalk are listed in TU>le 1. At asimuthal angle of 
0 and 90 degree the partide motbn is linear polarized. The part4cle 
motion with source at arbitrary aaimuthal angle did not show any 
particular p^ttoma. The partida motion pattern Is depend on the 
relative phase of the Green S function G„ and This phase is 
depend on the source receiver distance, degree of anisotropy and 
frequency. Ftom the particle motioa analysis we only see the initial 
partide motion b directed ateng the faat shoar wave direct witb 
large source receiver separation. But doe to the near field nature 
Of acoustic iog^ng this initial particle nsotion is hiud to obtain ac- 
coratdy for any practicle use. In order to answer the fioaj <)ue6tion 
we h«ve to seek another way to do the inverse decomposition. 

Ztivme Decomposition of Dipole Log 

In the field logging dtuations, we don't know the pHndple direction 
of the formation. The inverse problem b how to find the principle 
direction from the meaattpMnent at ar^trary admuthal angle. We 
assume two pair of Inline aad croGsline dipole moaAurtincnts are 
collected in perpendicular directions. The measured data is repre^ 
eented asadata vector V = (t7(r,), ufi,)^ The two inline 
dipole roeasuremenLs are Vj and vj***. The two crossUne dipole 
measucemenes are Vj, and y*wo. The problem beeomca how to 

find oximutKoI ftogl^ Mq<l inline meAaurtimentu aJojig tnc X micI 
Che Y ads from these four data veetera. By using Equation (5) and 
(6) we have: 



Med ftoin<7132388ll08>at 7/2210312:39:23 PMIEastemDayflgMM^ 



Sent by: CONLEY ROSE, P.O. 



713 236 aOOB; 



07/22/03 10:30; JfiiEBUL_#648i Page 4/5 



1 



Oecompotition of dlpole tog 



T 



(Vs-V,)dii^c<atf 

Equition (U) U not indep«a<kot oae bocOaue = - Vj^**- So 
it left ua witb three «qu4tiOn:» with three uqkoowra. Th« three un* 
knowia ^ aoglc «, V^. With a few algebm maiUpaUtions 

of £<)iiAtio& (11), (12) m4 (IS) we obtain: 



(11) 

(IS) 



V^(cantf-cot0) = 
V, « 



(15) 
(16) 

(n) 



notice that Equation (15) contalites d equ&tioA«, wKef« a b th« 
(iu(nbtt of aiea$ut«meata. So angle 6 is overdttermmed. We seek 
the least-^oare solution for Thb solution U given belo^i^: 



taad -1 = 0 



(18) 



0 be e^ay to obl&io by fiqdliig.the root of accoad-order poly do- 
m\M. Orud is known, V« tutd Vy can be cakaUted dir«C!ly frodn 
BqnAtlon (16) Md (17). So Inveise deoompooilioa U £tv>ed the 
anaJyticai expressions in the feast-sqMUe Seivse, 

Next we consider a poinericaJ ex&rti^lc. Austin chalk formation 
id coneidertid once agaiA. Flr^t we determine tbe rotsUed ft^gle 9 
by using Equuion (18). The inline and crOes^Iine dipob bg? ace 
obtained by u&ing Elqu&tton (5) and (6). Aad these dat» arc used 
to recover the rotated an^ Random itoisc uf added Into the inline 
and tbe cro^tine data to test 4eosittvity of the formula to nobc. 
The recovered ftn^les are lifted in IWe 2. If there is lo noise 
added into the data the rotated M^tcs ore recovered almo^ exaaly. 
With noifie the diifcrencefi between recovered aiig^o actu^ 
routed angle are les than 0.6 degree. When noise level increased 
to 10% the dUferenee^ arc less than US degree. Under the noise 
conditioQ the inverse decomposiUon formui^fitil) can ^ve very good 
eetiroation of the principle direcUon. Recovered waveforms dong 
the prindple direction are pbtted in Figure 5. The iaBne and 
cro&sline lagging daU at azimuthaJ angles of SO and 120 degree 
are used. The inGae and eroaarme waveform at aziroutba} angle 30 
can be found in Figure 2. The four traces shown In the plots are 
ori^nal data, recovered from notse free, S% noise and 10% noiae 
data. The waveforms along the principle directions are recovered 
Very well even under noisy condition. The inverse decompositioii 
fofiqulas are v«y compntatioanl efficient. It can be i»ed to pfoces» 
field acoustic lossing dat& in anisotfopic formation* 

Conclusions 

In this paper the analytic expressions are obtained for both for- 
ward and Inverse deoomposition of acoustic dipole logging data. 
These decottipoBittao foromlas-are valid because tbe principle of 
euperpositjon hold for linear wave propagation in anisotropic me- 
dia- These very ample analytic formuhvt can be directly applied 
to field logging data prooeiaiag. The ayoibetic acou»Uc dipOle logs 
ars used to demoe^ate these forward aod Inverse decompositLoni. 
The synthetics are calculated by 3-D finite difference method. The 
numerical examples arc shown that ipverse decomposition formula 
woHca vrry well with noisy data. The oarttcte motion of dlpob lo« 
in anisotropic formation is very complkatod. This is becaiiM the 
acflustle logging distance b In the near field. The puiicle motion 



with dipole source aligned with principle direciloft to lioeaj- polar- 
hal. When the dipole soiirce is at arbitrary asimuthal angle ibe 
Initial particle motion tend to point to the fast shear wave dliec- 
tbn. But this ioitial pol«v»ation can only have very limited use in 
the field dau situation. 
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Source direction 


HectUinearit 




(degree) 


shale (2ni) 


chalk (2m) 


chalk (4.5m] 


0 


1.00 


1.00 


1.00 


15 


0.9d 


0.97 


0.95 


30 


0.96 


0.93 


0.83 


45 


0.96 


0.93 


0.75 


60 




0^6 


0.B4 


T5 


0.99 


039 


0.95 


90 


1.00 


1.00 


1.00 



Table I: Polariaatian analysis of dipole log. The formations are 
Bakken shale and Aestin cbaik. The source-cecdvcr separuions 
are 2 sod 4.5 metens. 
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^cation 4ngle 


Recovered Mfcie 


(degree) 


noise (0%) 


noiiie (5%) 


Dol$e (10%) 


5 


5.00 


0.17 


5.S1 


15 


16.00 


13.<I3 


16.36 


30 


30.00 


30.44 


3] .22 




45^0 


4S.0d 


45.15 


60 


60.00 




59.24 


7S 


75.00 


74.82 


74-19 



Decomposition of dipole log 



tMt 2: Recovered nsution juig^e from IdUoa And cromljne data. 
ftaitdAm noise ^ added inlo the data. 




figure \: Cotti^Anson of the ftniu diflcruica synthetic with liitfiar 
decompoi&itian, Dipate source s at ayimatbiJ aagte of 45 decree. 
The Wfrvcfonn io the X direction il tthowQ. Tncb (a) finite dil&mice 
synthetic. Tku» (b) decompadtbn merdiag taEqaatba (1). Od 
left botrapic formatian. On right anisotro^c formation. 
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figure S; X-Y plane particle motion ofdipoie log in Austin chAlk 
formation. The 4^mutha| an^ of ihc dipole louree i& fa) Q (h\ 
IS. (c) 30, (d) 45. (e) 60. (f) 75. fe) flO, ^ ^ ^ 
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Figure 4: X'Y plane particle motion ofdipoie log in Austin chalk 
fonnaUoa. Th« «0UrC6- receiver diatarvoe id 4.5 meters. The aa- 
irauthal angle oT the dipole souica ia fa) 0. fb) 15. (c) 30 (d) 45 
(e) 60, (f) 7a. (g) flO. 



Figure 2i Projected dipole log In Austin chalk formation at different 
aaw|«thal angles by u^Ing Equation (5). On (eft inline djpole. On 
nght cro43i)ne dipole. 
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Figure 6; Recovered waveforme from 30 nnd 120 d^roe inline and 
cro»Dne dipde lo^ The formation is Ausein chalk. TVace (») 
original waveform along theX aioB. Traco (b) recovered from data 
without noise. Tfwe (c) recovered rim daU with 5% random 
noifie. Tf4ce (d) recovered from dat* with 10% random noise. On 
kfi the X direction. On rigiit the Y diraction. 
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